Nonlocality in kinetic roughening 
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We propose a phenomenological equation to describe kinetic roughening of a growing surface in 
presence of long range interactions. The roughness of the evolving surface depends on the long range 
feature, and several distinct scenarios of phase transitions are possible. Experimental implications 
are discussed. 
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"Suppose that we take a bin and gently and uniformly 
pour in granular material. As the material in the bin 
builds up we can identify a surface and ask the ques- 
tion, 'What is the magnitude of the fluctuation in the 
height of surface (measured from the base of the bin)?' 
Also of interest is the length scale of the surface fluctua- 
tions and how they behave dynamically as more material 
is added." [EJ And thus was born the Edwards- Wilkinson 
(EW) model for surface growth - a solvable linear model 
at the heart of our current understanding of numerous 
growth processes. A relevant nonlinear term, added to 
this by Kardar-Parisi-Zhang (KPZ) §-[§, brought to 
light the nuances of growth phenomena to the extent 
that the KPZ equation very soon became a paradigm, 
in particular for dynamic phase transitions. The appli- 
cability of the KPZ equation seems to encompass length 
scales from atomic level to macroscopic phenomena of 
every day life, but still a specter is haunting the field : 
why is the KPZ behavior not observed ]^)? 

Many of the experimental situations however involve 
complex processes which beg to go beyond the idealiza- 
tion as pouring of noninteracting particles. This is espe- 
cially true if medium or fluctuation induced interactions 
interfere with the process as for example in the several 
recently studied systems involving proteins, colloids or 
latex particles [^|-^), or in sedimentation. The major in- 
teraction one has to reckon with, as detailed numerical 
computations suggest P,[L0|, is the long ranged hydro- 
dynamic interaction. Are such long range interactions 
relevant for the roughness of the surface? This question, 
the absence of a formalism to handle such interactions 
in the growth process, and the elusiveness of the KPZ 
behavior, led us to propose a simple phenomenological 
model by focusing on the long range nature of the extra 
force. 

We like to develop a Langevin equation type descrip- 
tion where long range aspects can be simulated by a force 
at each point of the growing surface exerted by the par- 
ticles away from it - a hint to go beyond a strict local 
description. In the linear EW model, the growth is along 



the global normal to the surface without any overhang. 
The height h(r, t) at point r and time t satisfies the diffu- 
sion equation with an additional noise term. If, instead 
of the global, the local normal is favored, the KPZ (V/i) 2 
term is needed 0. This nonlinear term describes the 
lateral growth at a point as can be seen from the height 
profile • We now extend this physical interpretation 
and take the gradient (or its magnitude) as a measure of 
the local density of deposited particles. The long range 
effect is now incorporated by coupling these gradients at 
two different points. Based on this intuitive picture, the 
equation we propose is the following: 

dh{r,t) 



dt 



nV z h(r,t) +rj(r,t) + 



i J dr' i?(r')V/i(r + r', t) • V/i(r - r', t), 



(1) 



where k is the diffusion constant for the particles on the 
surface, and 7/ is a random space time dependent white 
noise of zero mean and (ry(r, t)r](r' , t 1 )) — 2A<5(r — r')5(t — 
t'). The kernel ?9(r) is of long range and, in principle, 
connected to the underlying interactions Ell ]. So, we 
take i?(r) to have a short range part Ao<5(r) and a long 
range part ~ r p ~ d , or more precisely, in Fourier space, 
tf(k) = A + X p k-P. Eq. @ then smoothly [§J| goes over 
to the KPZ equation for A p = 0. We show that this 
leading term introduced is sufficient to yield a new fixed 
point with continuously varying exponents, and different 
phase transitions not found in the KPZ problem. The 
connection with experiments is discussed near the end of 
this paper. 

A central quantity of interest in growth problems is the 
scaling behavior of fluctuation of the height, (| h(r, t) — 
h(0,Q) | 2 ) which on a large length and time scale has a 
scaling form | r | 2x T{\ t |/| r \ z ). Here \ is t ne roughness 
exponent of the growing surface and z is the dynamic 
exponent. These two exponents define the universality 
classes of roughening. 

At d — 1, for the local growth (i.e. KPZ) equation 
a disorder dominated rough phase is found for all Ao by 
several exact treatments providing % = 1/2 and 
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z = 3/2. The nonlinearity is marginally relevant at d = 2 
and for d > 2 there is a phase transition from a strong 
disorder dominated phase (x + z = 2 for all d) to a weak 
coupling phase where nonlinearity is irrelevant, i.e. a flat 
phase with 2 = 2. The perturbation theory is inadequate 
for the strong coupling phase at d > 2 due to the lack 
of a perturbative fixed point Jl5[ . Numerical simulations 
p6[ predict z sa 1.6 at c? = 2. The phase transition is, 
however, under control, with z = 2 Vd > 2 ]l7| , ^8[ , with 
a rather complicated critical behavior. ]l9| , |20| . 

A simple scaling analysis indicates that both A p and Ao 
are relevant for d < 2 at the Gaussian fixed point (EW) 
where one expects, \ = (2 — d)/2 and 2 = 2. This follows 
from the scale invariance of Eq.([j]) under the transfor- 
mation r — > 6r, i — > /i — > 6 x /i, whence k — > b z ~ 2 n, 
A -> b z - d - 2 *A, A -► &x+^-2 Ao and Ap _^ 6 *+x+p-2 A/)< 

Also for any nonzero A p with p > the local KPZ the- 
ory (A p = 0, and x + z = 2) is "unstable" under renor- 
malization and a non-KPZ behavior is expected. For 
2 < d < 2 + 2p only A p is relevant at the EW fixed 
point. In the following we adopt a dynamic renormal- 
ization group (RG )procedure. Our results show a new 
stable fixed point at d = 1, for any p > 0. Another inter- 
esting consequence of this nonlocality is the possibility 
of a stable fixed point at d = 2 for a certain range of 
p. The marginal relevance of nonlinearity in the original 
KPZ theory is destroyed. 

The renormalization procedure is most succinctly de- 
scribed through the Fourier modes momentum q and fre- 
quency to in terms of which Eq. (|l|) becomes 

h(q,w) = G (q 1 w)fo(q > u;) - (l/2)(2 7 r)- d - 1 x 
J dq'<L/i?(2q') q+ q h{q+, w+)/i(q_, w_)], (2) 

where, symbolically, AT± = AT/2 ± A"' with X = q or to. 
Here Go(q, w) = 1/(k<z 2 — iw) represents the bare prop- 
agator or the Green function for the diffusion equation. 
We follow the usual iterative perturbation scheme where 
h in the RHS of Eq. (|) is replaced by Eq. (|) itself 
upto 0{d 2 ). A convenient diagrammatic representation 
can be set up from this scheme and the renormalization 
of the various parameters can be obtained from appropri- 
ate vertex functions. We skip the details as they are very 
similar to Ref. pl| . In the subsequent renormalization 
procedure, we integrate out small length scale fluctua- 
tions over a momentum shell Ae~ l < q' < A to obtain 
the effective parameters for a similar equation but with 
a smaller cutoff Ae~ l , where A (set to 1) is related to the 
microscopic cutoff. A subsequent rescaling then restores 
the cut off to A. 

The effective propagator G(q, u) = h(q, w)/rj(q, u), 
gives the renormalization of tension k. The ef- 
fective noise, obtained from (h*(q,u>)h(q,ui)) = 
2AG(q, lo)G{— q, — u>), gives the renormalization of the 
disorder. Next we look for the terms contributing to the 



effective nonlinearity. Note that the RG transformation, 
being analytic in nature, cannot generate a singular term 
to renormalize A p for p > — 2. In fact there is no renor- 
malization of Ao either. A contribution to Ao could come 
from terms of 0(A$ 3 ), and a straightforward calculation 
pl[ shows that such terms do cancel each other p2] |. 

Following the above procedure we arrive at the flow 
equations for k and A as 



dn 

~dl z 

dA 



2- 2 



2) + 3/(1) 
Ad 



(z-d~2 x )A + 



A 2 K d 
4k 3 



$(2) 2 



(3) 
(4) 



where f(a) = dln{)(k)/d\nk \k= a , the (effective) expo- 
nent of -d(k) and Kd = Sd/{2-K) d , Sd being the surface 
area of a d— dimensional unit sphere. The flow equations 
for Ao and A p , having contribution only from the rescal- 
ing, are dX x /dl — (x + z — 2 + x)A x , {x — or p). The 
two parameters \ an, i z are chosen to keep k and one of 
X x invariant. 

In terms of U 2 = A\ 2 K d /K?, and R = U /U p , with 
the choice x + z — 2 (or 2 — p) and z equal to the ex- 
pression inside the big round bracket of Eq. [| the flow 
equations can be combined into two as 



dUp 
dl 



(2-d) 



U 



2d- 3. 



-rr 



U U P 



coUo + dUp], (5) 



and dR/dl = -pR, where c = {5d-6)(l + 2-P) -2d-9p, 
and ci = [(3 + 2- p )d - 6 - 9p]2" p . The equation for R 
rules out the existence of any off axis fixed point in the 
Uq and Up parameter space ( except for p = 0, when 
there is a trivial marginal fixed line.) 

There are only two sets of axial fixed points, SR= 
{Uf = 2d(d- 2)/ {2d- 3),uf = 0}, with x + z = 2, 

and LR = {U^ = 0, uf = Ad(d - 2 - 2p)/c x } with 
X + z = 2 — p. The first set (SR), with A p = 0, corre- 
sponds to the known KPZ fixed point, whose properties 
have already been mentioned. However we see a relevant 
perturbation U p which grows at this fixed point. The 
stable fixed point for d < 2 + 2p, with p > is LR, ex- 
cept for the region bounded by d = (9p + 6)/(2~ p + 3) 
and d = 2 + 2p. This excluded region is, like the KPZ 
case, an artifact of one loop renormalization p|. At the 
new fixed point LR 



2 + $, and x 



-P- 



(6) 



where $ = (d- 2- 2p)(d- 2- 3p)/[d(2~P + 3) - 6 - 9p)} 
This fixed point admits z < 1 (not unexpected for long 
range cases) but by virtue of the relation x + z — 2 — p, 
X need not be greater than 1, a requirement for ignoring 
higher order terms in Eq. [l]. At d = 2, the marginal 
relevance of Uo of the KPZ theory is lost and there is a 
stable fixed point (LR) for p > 0.194. 
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(A) d<2 
LR rough 



LR rough 



(B)2 + 2p>d>2 Ap 
LR rough 



(C)d > 2 + 2p 
LR rough 



LR rough 



LR rough 



(D)d<2 + 2p |(E)2 + 2p<d<2 1 T ) J > 2 



SR rough 



SR rough 



SR rough 



SR rough 



SR rough 



SR rough 



FIG. 1. A p vs Ao phase diagram. (A)-(C) correspond to 
p > 0, while (D)-(F) to p < 0. Thick line along the y-axis 
represents the LR phase, while the medium thick line along 
the x-axis SR phase. In (D) and (E) the phase is SR (KPZ) 
type for all A p 7^ 0. The dashed line in (B) and (E) represents 
a smooth phase, which extends over a region in (C) and (F). 

To discuss the surface morphology and the phase tran- 
sitions (23), we consider different values of d and p > — 2. 
See fig. 1. Also note that the invariance of Eq. [l] un- 
der h — > —h and A — > —A, is respected by the fixed 
point equations. Since the nonlinear term is like a force, 
the change in sign of A corresponds to a "push" -"pull" 
change or a growing to a receding surface case. We there- 
fore consider both positive and negative values of A p , and 
without any loss in generality, take Ao > 0. 

Case I: d < min(2, 2 + 2p): For p > 0, the stable fixed 
point, if exists, is LR with the dynamic exponent given 
by Eq. ^. Even if it does not exist in this one loop 
approach, still, from the flow, the phase is the "strong" 
disorder type. We call this an LR phase to distinguish 
it from the SR or KPZ phase. It is possible to have a 
transition between two identical LR rough phases (push- 
pull). The critical behavior is EW type if there is strictly 
no short range nonlinearity else it is KPZ (SR) type. See 
Fig 1A. In contrast, for p < (Fig ID), the LR is irrel- 
evant and the surface behavior is always SR (KPZ) type 
except for Ao = 0, when it is an LR phase. There is no 
phase transition for Ao =/= 0. 

Case II: min(2,2 + 2p) < d < max(2,2 + 2p): The 
phases are LR or SR depending on the sign of p (Fig 
IB, IE). For p > 0, the critical behavior depends on the 
strength of the SR nonlinearity, Ao- For small Ao < Ao c , 
the critical surface is a smooth one while for Ao > Ao c it 
is KPZ. There is no transition if p < and Ao ^ 0. How- 
ever, for Ao = 0, there is a LR rough to smooth transition 
for p < 0. See Fig. IE. 

Case III: d > max(2, 2 + 2p): For p > 0, the LR fixed- 
point is unstable. A small nonlinearity dies down yielding 
a smooth surface while large nonlinearity will produce an 
LR rough phase. Unfortunately the absence of a fixed 
point forbids any prediction of the behavior of the LR 



phase. The unstable LR fixed-point controls the tran- 
sition between the rough and the smooth surface with a 
dynamic exponent z c — 2+cpe+0(e 2 ) where e = d—2—2p 
and c = -2^/[(2 + 2p)(l + 3 2?) - 3 2^(2 + 3p)]. This is 
in striking contrast with the believed to be exact result 
of z c — 2 for the KPZ case. The phase diagram is shown 
in Fig 1C. For p < 0, there is a phase transition between 
a SR rough and a smooth phase only if Ao is less than a 
critical value as shown in Fig. IF. 

Experiments on colloids || have yielded a value of 
X = -71 which is also the value obtained from paper 
burning experiments [p4[ . These are taken as the ex- 
ponent for a driven surface (line in the d = 1 example). 
For the colloid problem, hydrodynamic interaction (HI) 
is important while in the paper burning experiment, it 
is possible to have a long range interaction through the 
microstructure of the paper. With this \ Eq. ^ at the LR 
fixed point in d = 1 gives p = — .12. At this point, it is 
difficult to conclude if this is the transient exponent seen, 
eventually going over to the KPZ value on large length 
and time scales (Fig. ID), or it is a true Ao = case. 
Other cases where HI is known to play a role, namely 
the deposition of latex particles or proteins, the experi- 
ments have not been done for roughness of the growing 
surface. We believe such experiments will shed new light 
on growth phenomena. 

For the KPZ problem, it is known that anisotropy of 
the substrate can lead to an overall irrelevance of the non- 
linearity in two dimensions [p5| . To see if anisotropy can 
have a major effect in the long range case, we now con- 
sider a variation of the problem where the long range in- 
teraction has different amplitudes in different directions. 
Restricting ourselves to d = 2, we take 

dh(r t) 

y = K\\df\h(r, t) + K±dlh(r, t) + t)(t, t) + 

I dr' ^ss,(r')ds!,h(r + r',t)d^h{r-r',t). (7) 

*=II.-L 

as the anisotropic version of Eq. [j]. In the isotropic 
case, r\ = i?||(r)/i?_i_(r) = 1 and r K = K\\/k± = 1 repro- 
duce eqn [j]. For simplicity let us concentrate only on a 
case of anisotropy in the long range part, with Ao = 0, 
Aj_(q) = \± p q~ p and A||(q) = A|| p <7~ p . An anisotropic 
scaling of the surface x± — > e l x± and at|| — ► e x\\ lead 
to r\ — > e 2 ( 1 ~^ l r\. For nonzero r\ the scale invariance 
consequently restricts Q = 1. For r\ = 0, this constraint 
cannot be imposed and the analytical tractability is lost. 
To avoid this complexity here, we take r x ^ and C = 1- 
The RG procedure follows as before, only a new flow 
equation for r\ is required. The recursion relations are 



dn± 
~df 



(z - 2)k ± + 



I6rl /2 2P 



r x \ 3pg±K±(r K +r\ + 2y/r K ) 



4rl /2 2P 



(1 



(8) 
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dr K 
dl 



1/2 

16 x 2P 



dg± 

~df = 2P91 



g 2 ± 2-P 



3pg±2~ 



16r, 



1/2 



(3r 



3r? 



4r,T(l 



2r A r K ) 



2r 3 / 2 



2ri 



i 2 i+P 



(9) 
(10) 



where g± = \\AK2Kj_ ■ For r x > 0, Eq. § has a 
fixed point with r* > 0, which is a continuation of the 
isotropic fixed point r\ = r K = 1 for p = 0. The im- 
portant fixed point for us is r* = — r\, which is phys- 
ical, from the stability requirement of the surface, only 
if r\ is negative. We consider only this anisotropic case 
here. The flow equation for g± now allows a fixed point, 
unlike the isotropic case discussed earlier. For small p, 
the anisotropic fixed point is at g\ ~ 8py/\ r\ |, with 
2 = 2 - p/2 + 0(p 2 ) from Eq. |. 

The effect of different signs of A p is to have opposing 
(push/pull) effects in the two orthogonal directions. In 
the KPZ case, they cancel each other producing a EW 
surface |25|]. In the LR case, we predict a new type of 
rough surface with z < 2 for p > though on the whole 
it may be flat but singular with x < 0. Surprisingly, this 
case seems to be better controlled in the RG approach 
than the isotropic case. This, in turn, calls for further 
studies of the r\ = situation entailing anisotropic scal- 
ing of space (C 7^ 1). 

In summary, we have proposed a simple phenomeno- 
logical model, Eq. |l| that incorporates, as a minimal 
model, long range interactions in growth problems. We 
have shown that any interaction decaying slower than 
\/r d makes the KPZ or the short range nonlinear case 
unstable, and asymptotically the surface will have differ- 
ent roughness with exponents depending on the power 
law of the interaction. The critical behavior in going 
from a growing to a receding surface can be of various 
types depending on the dimensionality and the strength 
of the interaction, as shown in Fig 1. Power laws de- 
caying faster than l/r d are suppressed by any local or 
short range nonlinearity yielding a KPZ like roughness, 
but when alone, it can produce a still rougher surface. 
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